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1 Introduction 



Perturbations of conformal field theories by relevant operators have been intensively studied starting 
with the work of Zamolodchikov [TJ [2] on integrable perturbations of conformal field theories. 
Numerous examples have been considered, but there are also a number of structural results, in 
particular the c-theorem of Zamolodchikov [3] that states that the central charge c cannot increase 
along renormalisation group flows, as well as the analogous g-theorem [3] for the boundary entropy 
[5]. Perturbations of conformal field theories also play an important role in string theory, for 
example for time dependent backgrounds, see e.g. [Q[7]. 

In the context of string theory also marginal perturbations are of significance. Most string 
theories of interest possess moduli, i. e. free parameters such as the size and shape of the background 
or the position of some D-brane, and these correspond to marginal operators in the two-dimensional 
world-sheet theory. Exact conformal field theory solutions are often only available at special points 
(in particular the rational points) in moduli space, and it is important to learn to control the theory 
away from these special points, i.e. after perturbations by marginal operators. 

Usually one thinks of the moduli as corresponding to exactly marginal operators, and then 
the renormalisation group analysis is, by definition, trivial. However, in the context of bulk and 
boundary perturbations, the situation can be more subtle. In particular, exactly marginal bulk 
operators (describing moduli of the closed string background) can cease to be exactly marginal in 
the presence of a boundary. If this is the case, they induce a non-trivial renormalisation group flow 
on the boundary [8]. 

In the analysis of perturbations by relevant operators, a first order analysis is usually sufficient 
(see e.g. [HI [TU] ) . However, in the context of marginal perturbations, it is often necessary to go to 
higher order in perturbation theory. The simplest example of such a situation is a single Neumann 
brane on a circle, for which the conformal dimension of the momentum eigenstates depends on the 
radius modulus. From the point of view of the world-sheet, the change in conformal dimension for 
this boundary operator does not arise at first order, but only appears at next-to-leading order in 
perturbation theory. 

Higher order conformal perturbation theory also plays a role in proofs of integrability of par- 
ticular bulk and / or boundary perturbations fTTl [T2"] . Conformal perturbation theory at higher 
orders was studied for particular models in [131 H US US]; general aspects of the pure bulk case 
were also discussed in [T7J, HE] • 

In this paper we make an attempt at a systematic analysis of conformal perturbation theory 
beyond the leading order. We begin by analysing which RG coefficients are scheme-independent 
(or universal) and thus can have a physical interpretation. (In particular, we show that this is the 
case for the coefficient describing the change in conformal dimension of the momentum fields on 
the Neumann brane.) We then outline a specific scheme — the position space minimal subtraction 
scheme — in which higher order RG coefficients can be calculated. This allows us to prove that 
the combined bulk-boundary perturbation problem is renormalisable at the quadratic order. While 
the minimal subtraction scheme is conceptually clean, explicit calculations of the RG coefficients 
are often rather cumbersome. We therefore also consider another, Wilsonian type scheme, to 
which we refer to as the 'OPE scheme' since the first nontrivial terms in the beta functions are 
given by OPE coefficients. This scheme has some conceptual shortcomings at higher orders but is 
computationally somewhat simpler. For the universal coefficients we are interested in, the result 
is independent of which of the two schemes we use (as we also verify explicitly). We can therefore 
determine the coefficients of interest (in particular the formula for the shift in conformal dimension 
for the momentum fields on the Neumann brane) in the Wilsonian approach. In the resulting 
formulae the universal quantities are expressed as integrals over certain correlation functions. As 
an illustration we also apply the formulae to an intersecting brane model on a torus, and again 
reproduce the geometric result. 
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The paper is organised as follows. In section 2 we discuss which RG coefficients in the boundary 
beta function are universal in the presence of marginal bulk perturbations (section 2.1). We then 
introduce the minimal subtraction scheme, both for pure boundary perturbations (section 2.2.1), 
as well as for the combined bulk-boundary problem for which we prove renormalisability at the 
quadratic order (section 2.2.2). We also introduce the Wilsonian scheme and discuss its advantages 
and shortcomings (section 2.3). Finally, we explain how the discussion can be generalised to include 
boundary changing operators (section 2.4). In section 3 these ideas are applied to two examples, 
the single Neumann brane on a circle (section 3.1), as well as a configuration of two intersecting 
Dl-branes on a 2-torus (section 3.2). Finally, we discuss in section 4 how our techniques for the 
calculation of higher order RG coefficients can also be applied to pure bulk or pure boundary 
perturbation theory. 

2 Bulk- boundary perturbations of BCFTs 

Let us start with a general discussion of the effect of perturbations by marginal bulk operators on 
boundary degrees of freedom. Generically such a deformation will induce a renormalisation group 
(RG) flow on the space of boundary conditions [8]. Under a certain condition, to be formulated 
precisely below, the induced boundary deformation is however scale independent to the first order 
in the bulk deformation parameter. In this case one can study how the set of boundary scaling 
dimensions changes with the bulk deformation. We will demonstrate that (despite no occurrence of 
RG flows) the RG technique is very useful in addressing this question. In particular we will derive, 
using the RG methods, general expressions for a first order change in dimensions of boundary 
operators along a bulk deformation. 

2.1 Universal terms in marginal bulk perturbations 

Before we give a detailed discussion we need to introduce some notation. Consider a boundary 
conformal field theory (BCFT) defined on the upper half plane 1H + = {(x,y)\y > 0} with complex 
coordinate z = x + iy. Let 4>k(z, z) be bulk primary fields with conformal weights (/i&, hk) so that 
their scaling dimensions are = 2hk- For a single (fundamental) conformal boundary condition 
we denote the boundary primaries by ip p (x), and their scaling dimensions by h„. Later we will 
generalise our discussion to superpositions of conformal boundary conditions. We will assume that 
the two point functions are normalised as 

5' 5 

{(l)i{z,Z)(j)j{w,w)) = , (lp p (x)tp q (y)) = ~~~To7T • (2.1) 

\z — w\ zn - 1 \z — W\ An P 

In particular, this means that we assume all fields to be self-conjugate; this is obviously not a real 
restriction, and our analysis can easily be generalised. The operator product expansion (OPE) for 
pairs of bulk and boundary operators has the form 

U^z^jiw,™) = J2 C ^ k \ z - w \ Ak ~ Al ~ A] Mw,w) + ... , (2.2) 

k 

Mx)4> q (v) = E Vfo - x^-^-^Mv) + ■■■ (v > ■ ( 2 -3) 

r 

Finally, when a bulk operator approaches the boundary it can be expanded using the bulk to 
boundary OPE 

<j) k (x + iy,x- iy) = Y, B k P(2y) h ?- A *ip p + ... . (2.4) 
p 
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With these preparations, let us now consider a perturbation of the given BCFT generated by 
the Euclidean action perturbation 



x ' ' A >'- 2 ^ f fdxdy4> k {x iy )+Y,l hp ~ 1 ^ fdx^ p {x) . (2.5) 
k ^ p J 

Here X k , fi p are the dimensionless coupling constants of the respective operators, and I is a renor- 
malisation distance scale. Up to second order in the coupling constants, the beta functions have 
the following general form 

P k = Vk\ k ■ >>;;A'A' • ... (2.6) 

ij 

i qr iq 

Here, yf. = 2 — and y q = 1 — h q are the bulk and boundary anomalous dimensions, respectively. 
The omitted terms stand for higher orders in the coupling constants. A general property of any 
local RG scheme is that the bulk beta functions are independent of the boundary couplings. 

It has been known for quite some time |19j (see also [20j) that in a particular renormalisation 
scheme the coefficients for bulk theories are given byC& = tt(7-A where Cm* are the bulk OPE 
coefficients from (|2.2|) . The same scheme, to be discussed in more detail in section 2.3, can be 
easily adapted for theories on the half plane. The coefficients T> p r coincide then with the boundary 
structure constants D qr p (see e.g. and for the coefficients B\ we have B\ = \Bp , where Bp 

are the bulk to boundary OPE coefficients from (12. 4p . see [8]. 

Consider now the case where we perturb the BCFT by a single bulk field (f>(x, y) with a coupling 
constant A. Furthermore we want to assume that the bulk beta function /3^(A) vanishes. In this 
case, even in the absence of an initial boundary perturbation /U^ are = 0, a boundary renormalisation 
group flow can be triggered by the terms B^\ in the boundary beta function. Such boundary terms, 
however, are in general not universal. For example, if the induced boundary fields are all relevant, 
i.e. y p > 0, then the corresponding terms in the boundary beta function can be removed by a 
coupling constant redefinition 

^ pp = + ^A . (2.8) 
Vp 

The above coupling constant redefinition looks peculiar in that jl p is not proportional to \i v . It has, 
however, a simple meaning. Let Z = (e is ) be the renormalised partition function^] of the perturbed 
theory (|2.5p . We have 

d\nZ \ (dlnZ\ v g fd\nZ\ (8\nZ\ 

) {ilP} V d\ ) {llP} ~ If y p { dp? J x 2^ yp { dfl p ) x > I ■) 

where the partial derivatives on the left hand side are taken with the boundary constants fi p or fl p 
held fixed. The identities (|2.9p mean that, after the redefinition (|2.8p . the bulk coupling constant 
A couples to a re-defined field 

B p 

4>(x, y) = 4,(x, y)-J2~ ■ ( 2 - 10 ) 



1 Here we are talking about the partition function on a disc. The passage from the half plane to the disc is 
straightforward because the bulk theory is conformal. 
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Part of the renormalisation procedure amounts to defining the operator coupling to A on the half 
plane so that the correlation functions involving that operator are distributions (and thus integrable 
in any bounded region on the half plane). In the interior of the half-plane the resulting operator 
must coincide with the bulk operator (f>(x,y), but in general extra subtractions may be required 
at the boundary. The redefinition (I2.10p stemming from the change of scheme (12. 8p reflects the 
natural ambiguity in defining such a fully subtracted operator extending <p(x,y). 

Suppose now that all terms linear in A can be removed in this manner. Then the resulting 
boundary beta functions have the form 

^ = J2D p q (X)fi q + 0(fi 2 ) , (2.11) 
<? 

where 

D p q (\) = y P 5 p q + A £* q + 0(A 2 ) with £* q = £? q - 2 £ V* qr) 5 ( 2 .12) 

and 23?.) = \^D\\ + T^rq)- Now that the boundary beta functions are all proportional to the 
boundary coupling constants we can treat the boundary perturbations (at least those of the relevant 
operators) infinitesimally to read off the dimensions of the boundary operators in the deformed 
theory. More specifically, we claim that the eigenvalues of the matrix Dq(\) are to be identified 
with y = 1 — h, where h is the scaling dimension of the boundary operator in the deformed theory, 
and the corresponding boundary primaries are the eigenvectors of Dq(X). To leading order in A, 
the matrix Dq(\) can be diagonalised by the transformation 

( £ p 

^-E^ + A/P)^, where f» = rf? for P / 9 (213) 



for p 



The corresponding primary fields ar^| 



V; P [A] = ^ - AV — ^—ipq , (2.14) 



and their anomalous dimensions are 



y P [A]=y P + A^ p . (2.15) 



We further claim that the quantity specifying the dimension shifts 

Vr 

is scheme independent. To see this we consider a coupling constant redefinition of the form 

^^^^^ + ^^7 + ^^ + ... . (2.17) 

qr q 

Under this redefinition the coefficients in the beta functions (12, change as 
V P s ^V P s + d p rs) (y r + y s -y p ) , 

£jr " 8 lr ~ 2 E VP r* bS ~ 2 E ^rsfiVr + V. ~ Vp) + 2 £ ^BJ + e^r ~ V P ) • (2.18) 

s s s 

It is straightforward to check that under the transformations (|2.18|) the quantity (|2.16|) is indeed 
invariant. 



The fields ip P [A] are defined up to adding a multiple of Aip p . 
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2.2 Computation in a minimal subtraction scheme 



In the following we want to explain in detail how these coefficients — in particular (I2.16P — can 
be calculated explicitly. We shall first study this question in a minimal subtraction scheme. In 
order to make sense of the formal perturbation series we shall use a point splitting regularisation. 
In particular, we require that any two perturbing bulk or boundary fields do not approach each 
other closer than a cut-off e, and that the perturbing bulk fields only approach the boundary up 
to a distance e/2. Before specialising to the bulk and boundary situation discussed in the previous 
section, let us first discuss some generalities of renormalisation. For brevity we consider only 
boundary perturbations but that is inessential for the points we want to make. 



2.2.1 Generalities of minimal subtraction schemes for boundary perturbations 

Let us consider a perturbed BCFT action 

S = Sbcft + XX l dx M x ) > ( 2 - 19 ) 
v J 

where (j^ are the bare coupling constants. Let I be an infrared distance scale at which we wish 
to renormalise the theory. In terms of the renormalised dimensionless coupling constants /i p , the 
same Lagrangian (|2,19|) can be expressed as 

S = Sbcft + Yl ryp » P [ d *M x ) + s a , (2-20) 
p J 

where y p are anomalous dimensions of the fields ip p (x), and 5 c t is a counterterm action. Perturbation 
theory generates integrals of the form 

...dxi. ..dx n ip Pl (xi)ip P2 (x 2 ) ■ . .ifj pn (x n )Y[8(\xi ~Xj\- e) 6(L - \xt - xj\) , (2.21) 

i<3 

where we have also introduced an infrared regulator L. The above expression is to be understood 
in the operator sense, i.e. inside a correlator with arbitrary other insertions. The product of fields 
in (|2.2ip can be expanded in terms of a complete set of local operators as 

ip P i{xi)i>p 2 {x2) •• -VwOcn) = ^C^_ jPn {xx,...x n -x)^ A {x n ) , (2.22) 

A 

where we have arbitrarily chosen the point of insertion on the right hand side to be x n . If the 
OPEs of the conformal families of the primaries ip p close on themselves, we can take for ^ a the 
fields tp p and their conformal descendants. In conformal field theory the expansion (|2,22p always 
converges [22] unlike in massive QFTs for which the OPE may be merely an asymptotic expansion. 
Substituting (|2.22p into (|2.2ip we obtain expressions of the form 

^C p A lv .. iPri (e,L) f dx^A(x) , (2.23) 
A J 

where 

C£ u ^ Pn {e,L) = J. . Jd Xl . . .dx^CpX^ixu . . .x n ^ 1 )f[9(\x i -x j \-e)9(L-\x l -x j \) . (2.24) 

i<j 

The integrals (|2.24p are finite because of the cut-offs. In the limit e — > the coefficients C pi 4>n (e, L) 
diverge with the divergences coming from the regions of integration in which two or more insertion 
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points xi, ...x n collide. In fact, if k operators ip Pl , . . . ,ip Pk come close together (with the other 
insertions bounded away from the point of coincidence) to produce an operator tyg, the leading 
divergence has the form 

Ce y n + '" +Vp k~ys ]jis+vp k+1 +-+vv n -VA ^ (2.25) 
where C is some numerical constant. Here we have assumed that the resonance condition 

Vpi + ■ ■ ■ + V Pk - VS = (2.26) 

does not hold; otherwise the corresponding divergence is logarithmic. The above reasoning follows 
essentially from dimensional counting, as well as from the locality of the OPE, ensuring the inde- 
pendence of the expansion (12.220 from L. If we only perturb by marginal or relevant fields, y p . > 0, 
then divergences can only occur if also ^> s is relevant, i.e. ys > 0. Assuming that the OPE is closed, 
^5 is then one of the perturbing relevant or marginal primary fields ip p . In the 'minimal subtraction 
scheme' we are using here, we only introduce counterterms for actually divergent contributions; the 
above reasoning then implies that the scheme closes on itself. 

The divergences arising when k operators come together first emerge at order n = k in perturba- 
tion theory. One expects that they can be canceled by local (L-independent) counterterms. These 
counterterms then also cancel the non-local sub divergences (|2.25p that appear at order m > k in 
perturbation theory. Thus we only need to deal with the case when n = k, in which case ^fs and 

must have a non-trivial two-point function, and hence ys = VA- Then the coefficient (|2,25p 
is independent of L, and hence converges when L — > 00. Note that the lower order counterterms 
may also contribute to the k = n divergence when the counterterm insertion from the order I < n 
comes close together with n — I fields ip Pi . The same dimensional reasoning however tells us that 
the final coefficient must again be independent of L, and the remaining divergence can be canceled 
by a local counterterm. 

The above discussion should however not be taken to be a recursive proof of renormalisability 
of conformal perturbation theory. One problem that needs to be tackled is the classical problem 
of overlapping divergences; in the case at hand this occurs when k points come together with a 
subset of I < k points coming together much faster than the remaining ones. The associativity of 
the OPE in conformal field theory should be the key property ensuring the consistency in dealing 
with overlapping divergences, but we have not attempted to work this out in detail. However, we 
will see in the concrete examples of the next subsections how the above discussion can be made 
more rigorous. In particular we will prove the renormalisability of conformal perturbation theory 
at the next-to-leading order using analytic properties of conformal blocks. 

In order to illustrate these ideas, let us now consider an integral that emerges at second order 
in perturbation theory, 

l^r^^'fci fdx 2 Mxi)M x 2)0(\x 1 -x 2 \-e)e(L-\x 1 -x 2 \) . (2.27) 
p,q 

The product of the two boundary fields can be expanded via the OPE (j2.3|) . Performing one of the 
integrals, we see that we get ultraviolet divergences of the form (in the limit L — > 00) 

where the summation runs only over those indices p, q, r for which y p + y q — y r < 0. In particular, 
y r > 0, and thus only relevant primary fields ip r contribute. As before, we have also assumed here 
that there are 'no resonances', i.e. that none of the expressions y p + y q — y r vanishes. Then only 
power divergences occur at this order. 
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(2) (2) 

The above divergences can be canceled by adding a minimal action counterterm S^t = ~^div- 
Equating the two expressions (12.190 and (I2.20p we obtain up to second order in the coupling 
constants 

1 D pq r (e\Vp+yq-y 



^ = r* M r + £ £ ^ T yr rw , (2-29) 



(2) 

where Ir is the set of pairs of indices (p, q) for which y p + y q — y r < 0. Differentiating both sides 
of (|2.29p with respect to / with fixed \f B we obtain 

l-^ r = f3 r (ri = yr» r ■ (2.30) 

Thus the beta functions are linear in fj,. It is easy to see that this property continues to hold also 
at higher order in perturbation theory, as long as the divergences are power-like. 

On the other hand, if we have a non-trivial resonance at lowest order, i.e. if y r = y p + y q , then 
formula (|2.28[) takes the form (we are assuming for simplicity that there are no other divergences 
at this order) 

S£ = -^ w r ln( £ /0r^>V P JdxMx) , (2-31) 
where we cut off the divergent integral in the infrared region at the renormalisation scale /. Intro- 

(2) (2) 

ducing a counterterm S^. = —S^ we then obtain a beta function for the coupling // 

P r = y rf i r + £> M W • (2-32) 

More generally, in the minimal subtraction scheme at hand, the nonlinear terms in the beta func- 
tions all come from resonances. However, in general not all resonant terms are universal. 

2.2.2 Minimal subtraction scheme for bulk-boundary perturbations 

After this interlude we now return to the case of interest, namely the description of the minimal 
subtraction scheme for bulk-boundary perturbations. In fact, the above discussion generalises in 
a straightforward manner to include an additional perturbation by a bulk field. For simplicity of 
presentation, we shall assume that the bulk field cp(z, z) is a spinless relevant or marginal primary 
field of scaling dimension A = 2 — with y^ > 0. As we shall explain below, the bulk-boundary 
perturbation at the next to leading order in perturbation theory is then renormalisable. We will 
specialise to the situation where the bulk field is marginal {y^ = 0) later. 

At the linear order in the bulk coupling A the divergences in perturbation theory only arise 
from singularities as the bulk field approaches the boundary. These are described by the bulk to 
boundary OPE. If there are no boundary fields for which h p = A — 1 we have power divergences of 
the form 

^ = -^^0)— (-3) 

where 

= {v\y P > ■ (2-34) 

In the minimal subtraction scheme the counterterm is then simply = — sQ\ If there is a 
boundary field for which B^, p ^ and the resonance condition y p = y^ is satisfied, we have a 
logarithmic divergence which results in a universal term linear in A in the boundary beta function 

m 

f3P = y pfJl P + X^f + ... . (2.35) 
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In the case when the bulk perturbation is marginal the resonance condition requires that the 
boundary field is also marginal. 

We shall, in the following, always assume that the resonance condition is not satisfied, i.e. that 
Up 7^ y<t>i this is for example true in the context of section [XT] where 2/0 = and y p > 0. Then the 
counterterm is given by S^- = — S^l. 

At the next order in perturbation theory we encounter the integral 

(2,6) 

where z = x + iy and R is an infrared regulator. The quantity in the square brackets in (|2.36l) 
can be expanded in local boundary fields as in (|2.22|) and (|2.23p : the coefficients of these fields 
can be expressed in terms of certain integrals (see below). By the same arguments as above, only 
coefficients of (primary) relevant fields can be divergent as we send e — > 0. 
More precisely, the coefficient with which the primary field ip p will appear in (|2.36p equals 

P q = JJdxdye(y-e/2)e(R 2 -(x-x') 2 -y 2 )(4>(z,z)4, q (0)^ p (w)) . (2.37) 

Using the Mobius symmetry, the correlation function appearing in this formula can be written as 

{4>(z, z)ip q (0)if) p (oo)) = \ z \vt-yp+v*- 2 n s+y *- 2 (1 - r? )( 2 -^- 5 )/ 2 Y(rj) , (2.38) 

where 

»7 = 1 — J , S = -(4-y p -y q - y<f> ) = -(h p + h q + A), (2.39) 

^ A B 4 , A DaA p e i ?toA-v++ 1 )F£r (r?) (Re*>0), 
Y(v) = < „, , HQP (2.40) 

^V^^-^W (Rez<0). 

Here the index A runs over all conformal primaries whose conformal families appear in the in- 
termediate channel. The conformal blocks (jl) have a branch cut along the real 77- axis from 
— oo to 1 and are normalised so that F A (i]) ~ <q h A-5 with coefficient 1 as T] — > 0. The conformal 
blocks entering the function Y{rj) are defined on opposite sides of the branch cut for Rez > and 
Re z < 0. The analyticity in z implies that the values of Y(rj) in the lower half 77-plane are obtained 
by the analytic continuation in a clockwise direction from the upper half-plane |23j . 

Passing to polar coordinates z = re 1 ^ and using (|2.38p we can rewrite (|2.37|) as 

71— 0* R 

p = [ d-d I drr v «- y ' +v +- 1 ri Hv *- 2 (l - rj^-^-^^Y (rj) , (2.41) 



and 



where 



tf, = arcsin (^) , r 2 ( V ) = ^ 1] . (2.42) 

Since 77 = 1 — e~ 2% ® depends only on d we can perform the integral over r. In the remaining integral 
it is convenient to change the integration variable i? to 77. Altogether we then obtain 

i f r / n — 1 \ ^pi/ 2 1 

P q = / d^e^flL—^j -RCn rf +v +- 2 (l-rj)- ( *++W*Y(jj) , (2.43) 

2 Cpq J L V V ) * 

C{e/R) 
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where 

C Pq = yq-y P + y<t> , (2-44) 

and the contour of /^-integration is a segment of the circle of radius 1 centered at r\ = 1 and oriented 
clockwise 

C{e/R) = {r, = l- e" 2M? , 0* < < tt - 0*} . (2.45) 

Obviously, this expression only makes sense if £ pg 7^ 0. In the resonance case (i.e. for Q pq = 0) we 
have instead 

(^)res = \ J drt In f ^- 2 (1 - ^"to+^yfa) . (2.46) 

C(e/«) 

The divergences of Iq and (if ) re s i n the limit e — ► can now be analysed using well-known 
properties of conformal blocks. There are two kinds of divergences that will be important to us: 
those that come from the region of integration rj ~ e/R — > where the bulk operator approaches 
the boundary far away from the point of insertion of tj)q\ and those that arise when the bulk field 
approaches the boundary in the vicinity of the boundary field ipq. As we have argued before (and 
as will become clear below) the former divergences are canceled by the contribution from the lower 
order counterterm S<£ , while the remaining divergences have the power e^ pq . To see this, we use 
the asymptotics — see (|2.40p 

for -> Y{r]) ~ ^ VV e ' !(M " ! '* +1 V" !/A " i + • • • , 

A 

for -> vr Y{r/) ~ ^ B/V^ ( "" M " 1) 1 1 " llA " i + • • • ( 2 - 47 ) 

in (|2.43p and ()2.46p . and then perform the rj integrals in the vicinity of rj = 0, i.e. from r] = ij^ up 
to some intermediate cut-off £. This leads to 

/P = + - V [ B^iDqAP + D Aq P) y^y AR y A+ y q -y p + Q/^-^+lO ( 2 _ 48 ) 

^\-2{y<i>-yA){yA + yq-y P ) J 



.4 



where C|", (Cq) res , and (/f )res are some constants independent of e and .R (that come from the 
evaluation of the primitive function at £, as well as from the remaining part of the integral). Since 
we are only interested in divergent contributions in e, we may restrict the fields A to be relevant 
primary fields in A E I^ 1 '. Furthermore, we can ignore all the subleading terms 0(e y <t>~ VA ^) since 
they vanish in the limit e -> oH 

Now we want to show that the divergent terms in the sum in (I2.48D are precisely canceled by 
the lower order counterterm given in (12.330 . At order Xfi q — recall that /i 9 is the coupling 
constant corresponding to ip q — the counterterm leads to the contribution 

Xn q r y "~^ Y ^"""oT-^Z — \ f fdxdx' 6{\x-x'\-e)6{R-\x-x'\)ip s {x)^ q {x'). (2.50) 

rm ^[y<t> ~ Vs) J J 

sEl W 



3 This is obvious for tja < 1. In a unitary BCFT tja — 1 corresponds always to the identity operator Q which does 
not have a descendant operator at level one since L-\Q — 0. 
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Again this can be expanded in terms of local boundary fields, and the divergence in the coefficient 
of ip p equals 



E 



B%(D q8 r + D aq P) 

L 2 (y<p - y s )(y s + y q - y P ) 



(2.51) 



Here we have assumed that there is no resonance among the boundary fields, i.e. that y s + y q ^ y p 
for any s € I^ 1 '; if there is a resonance, i.e. y s + y q = y p , then ()2.5ip has to be modified in the 
obvious manner. 

In either case, by comparison with (|2.48|) . it is now clear that the contribution (J^ ') q cancels 
precisely the divergent part of the sum in (|2.48|) . and similarly for the resonant case (, pq = 0. Thus 
the divergent contribution only comes from the first term in (|2.48|) and (|2.49[) 



1 q 



(Iq)rcs = {Iq )res + {{I c t )g)rcs ~ (Cg) res ln(e/Z) 

where C q and (C q ) Tes are coefficients that can be obtained by taking the limits 

(CP) res = \imed e (P q ) res . 



ct 



C p q e^ 



as e 



c\ 



lim e ed t P , 



Using the explicit expressions (|2.43p . (|2.46p and (|2.5ip we finally obtain 

% 



ci 



lim 



pq 



C(e/R) 

B/(D qs p + D sq p) {R\y°+y«-yp 

e 



(2.52) 



(2.53) 



(2.54) 



+ Jf?u 2 (y s + y q - y P )Cpg 



b^(d, 



qs P D sq P) 



+ s ^d 2 {ys-y<t>){y s + y q 



(Cf)ree 



lim 



I 

2 

C(e/R) 



E 



B^ s {D q P 



Dsq P ) 



2(y s - y<t>) 



jl\y s ~y<t> 

€ 



(2.55) 



It is worth noting that the contours of the ^-integrations in (|2.54[) and (|2.55p can be deformed 
provided the ends of the contour are held fixed and the cut is not crossed. In particular one can 
deform C{e/R) to run infinitesimally above the cut to r\ = 1, and then infinitesimally below back 
to the second endpoint near rj = 0. 

Given that the conformal blocks only have singularities at r\ = 1,0, oo with standard asymp- 
totics, the quantities C q , (C q ) Tes defined in (|2.54|) and (|2.55|) are finite and independent of R. This 
essentially provides a proof that at order XfjP the divergences can be canceled by local countert- 
erms which are linear combinations of relevant operators. Together with the well known results at 
orders A 2 and \jp\x q we have thus shown that a generic perturbation by relevant or marginal bulk 
and boundary fields is renormalisable at the quadratic order in the couplings^ It should also be 
possible to extend the analysis to higher orders in perturbation theory, but we have not attempted 
to do so. 

If the resonance condition Q pq = is not satisfied the divergence of I q is power like, and in the 
minimal subtraction scheme there are no terms of order \fi p in the beta function (5 V . On the other 
hand, when the resonance condition is satisfied for a pair (p, q) we have a universal term (c/. (|2,31|) 
and (12321) ) 



0P = y p ^-(CP) Ies X^ + ... . 



(2.56) 



4 It should be clear from our analysis that the different technical assumptions, namely that y s 7^ y$ for any 
boundary field ip s , and that we only have a single bulk field, are not crucial for the argument. 
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For a marginal bulk perturbation y^ = 0, the p = q case (and only that one) is always resonant so 
that we have in the notation of section 12.11 

{o (&p) ™ tr/ q . (2 - 57) 

Moreover, assuming as in section [27T1 that all y p > 0, we have in the minimal subtraction scheme 

2^ = 2^ = 0, B5 = 0, (2.58) 

and therefore by (|2.16j) 

% P = -(C p p U, (2.59) 

where (Cp) res is given by formula (|2,55p for = and p = q. As we proved in section 12.11 
the quantity Sq is scheme independent; we have therefore managed to obtain a description of this 
universal quantity in terms of conformal blocks — this is the main result of this subsection. 

2.3 Computation in a Wilsonian scheme 

It is instructive to compute (|2.16|) also in a different, Wilsonian type, renormalisation scheme - 
the one of [201121} [8]. We will refer to this scheme as the 'OPE scheme' for the reason that the first 
nontrivial terms in the beta functions are given by various OPE coefficients. This scheme is often 
employed in conformal perturbation theory at the leading order. One of the advantages of this 
scheme is that in the presence of a nearby infrared fixed point in theory space, the corresponding 
coordinates are nonsingular near that fixed point. Another attractive feature is that formulae for 
universal quantities, such as the dimension shift (|2.16p . can be obtained quite easily in contrast 
with the minimal subtraction scheme. On the other hand we will see at the end of this section that 
the scheme has some pitfalls when applied to computing non-universal quantities at higher order 
in perturbation theory. 

In the OPE scheme the theory is also regulated by a point-splitting cut-off e, just as in the 
minimal subtraction scheme of the previous section. It is however convenient to introduce the 
infrared regulator slightly differently: we introduce a cut-off whenever two coordinates X-^^ Xj £1X6 
separated by a distance larger than L, and whenever there is a bulk operator at a distance y > L. 
The dimensionless couplings are now introduced by using the UV cut-off scale itself, which is 
understood as a fundamental UV scale (lattice spacing, atomic or molecular scale). Thus we have 

SS = e Ak ~ 2 X k J I dxdy </> k (x, y) + e h "~ 1 pP j dx ^ p {x) . (2.60) 

k ' "' p 

Note that one can also include in (I2.60P irrelevant operators ipA with y^ < 0. Their contributions 
will be relatively suppressed as e~ VA but one may worry that in the perturbation expansion they 
will lead to contributions more singular than this suppression factor. We will see that, although 
there is no need to introduce irrelevant operators at the leading order, they are sometimes necessary 
to be taken into account at higher orders in perturbation theory. 

We will confine ourselves throughout this subsection to the case of a single marginal bulk field 
<j)(x,y) (A = 2). In this case the terms in the perturbation expansion we are interested in are 



e ss = 1 + A JJdxdy<t>(x,y)6(y - |) + ^V^/i p Jdx^ p (x) 



« V» I ldxdyt>(y- e -) jdx'4,(x, y )-ip r {x')e(L-\x-x'\)e{L-y) (2.61) 



P 

.^e-yp-Vi^P^l I ldx 1 dx 2 4> P (x 1 )'ip q (x2)9(x2-x 1 -e)6{L-\x 1 -X2\) + ... • 

pq 
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The cut-off variation ed e e ss can be computed assuming that the coupling constants depend on the 
cut-off via the couplings themselves according to 



(2.62) 



where (3 V are the beta functions (|2.7p . In the OPE scheme we now vary e ss with respect to e, 
i.e. we compute ed € e ss , and demand that the variation vanishes at the leading order in e. This 
reflects the main principle of the Wilsonian renormalisation group approach, namely that the 
renormalised quantities must be independent of the UV scale. The resulting equations fix order 
by order the coefficients of the beta functions. The linear terms in the beta functions are always 
scheme independent with the coefficients given by the anomalous dimensions. It is easy to check 
that at the linear order in jjP the equation 



ed e e ss ~ 

is satisfied automatically. The equation arising at the linear order in A fixes 



(2.63) 



(2.64) 



where B^ are the bulk-to-boundary OPE coefficients (12. 4h (see [8]). At the quadratic order in the 
boundary couplings one obtains the well known expression 

V? B = A-/ , (2.65) 
where D rs p are the boundary OPE coefficients (12. 3p . Finally, the equation at order \fi q is 







r e L ~ y i f f , e , 
o Xfi9 1 T~ 1 1 dxdx 2^ q ^ X " > 6<kL ~\ x ~ 



2 J J 

e ~y P ~y q 

+ 2^ B 4> P 1 1 dxidx 2 ip q (x 1 )ip p (x 2 ) 0(\xi -x 2 \ -e)9(L- \x\ - x 2 \) 



(2.66) 



The first line in the above expression came from applying ed 6 to the cut-off function 9(y — |) while 
the second line came from the lower order term (I2.64|) . The coefficients are formally obtained 
by taking correlation functions with the operator ip p inserted at infinity 



e dx(^(x,-)tj}q(0)tp p (oo)) 



£~ Vr B/ J dx (^(0)W(z)W(oc)} 9(\t 



(2.67) 



The above expression is formal because the limit may not exist. The integrals in the second line in 
(|2.67p can be evaluated explicitly 



I eW-y-* B/ f dx (^(0)V r (x)^(oo)} 0(\x\ - e) 



D (qrf 



(Vg -Vp + Vr 



L 



Vq-Vp + Vr 



ee H(e/L) . 



(2.68) 
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To study the convergence we rewrite expression (|2.67p via conformal blocks using ([2.380 



OPE 



lim 



d vv s -y^yp- 2 (l - r])- 5 ^(r] - l)(yp-v«)/ 2 Y(r]) - H(e/L) 



C'(e/L) 



where 



C'(e/R) = { V = 1- eT 2l \ < < < tt - <} , < 



1, (\-iel1L 

In ; — 

1 + ie/2L 



(2.69) 



(2.70) 



is a segment of a unit circle around r] = 1 oriented clockwise. If we now evaluate £ p p of (|2.16p in 
the OPE scheme, using (|2.69p for p = q, as well as (|2.65j) and (|2.64p . we obtain the same expression 
in terms of conformal blocks as given in (|2.55p and (|2.59p . Thus is indeed scheme-independent, 
as we have argued before. In terms of correlation functions, it can now be written as 



lim 



- /dx(^,-)^(o)^M)-^D 



(pr) 



Vr 



(2.71) 



2L 



and using the variable t] = 1 — e 2ie in the integral (|2.69p . we can also obtain 

4>p 

TV— 5 



Writing 5 

another, perhaps more elegant expression for £ p 



lim 

5^0 



t) 



Vr 



28 



(2.72) 



where the bulk field insertion runs over a semicircle of radius one around the boundary insertion 
ip p (0). Note that there is nothing special about the radius being one, since 



(0(e^)^(O)V P (oo)) = p 2 <<Kpe"Vp(0)W(oo)) 



(2.73) 



for any p > 0. 

Let us now come back to the expression (|2.69|) for p ^ q. Substituting the asymptotic expansion 
(|2.47p into (|2.69|) we find that although the most dangerous divergences, associated with the leading 
contributions of relevant primaries in (I2.47|) . cancel out, there may be divergences coming from 
terms in ()2.47p associated with irrelevant fields. More precisely there are additional divergences in 
(£^ 9 )ope from the region near r/ ~ whenever there is an irrelevant primary (/)a( z ) in the theory 
such that 

WO and {D qA P ^ or D A / ^ 0} and y A + y q > y p , (2.74) 
or whenever we have a relevant primary tp r (z) such that 



B/ ^ and {D qr p ^ or D rq p ^ 0} and 



y r + y q > y P + 1 • 



(2.75) 



In the last case the irrelevant field causing the divergence is a descendant of the primary ip r . From 
the point of view of the minimal subtraction scheme of the previous subsection, this problem does 
not arise since the conditions (|2.74p and (|2.75p imply that there are no divergences from the region 
where the bulk field approaches the boundary insertion. In fact, the additional e — > divergences 
in the OPE scheme come directly from the extra divergent factors of e included in the action (|2.60p . 

The situation can be mended if we include in the original perturbed action (I2.60P also irrelevant 
fields, and introduce their beta functions by requiring that ed e e ss ~ at the subleading orders 
in e. Then the function H{e/L) is modified accordingly to include more divergent terms that 
cancel out the divergences coming from the integral in (|2.69p . Although this resolution looks quite 
natural from the Wilsonian point of view, the whole scheme becomes quite unwieldy for practical 
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applications whenever (j 2 . T4 [) or (|2.75p happens. Note, however, that these extra divergences do not 
appear for universal quantities like E*L. Thus, as long as we are only interested in these quantities 
we can (and will) use the technically simpler OPE scheme. In particular, we will use this method 
to compute analogous quantities for pure bulk and pure boundary perturbations in section 4. 

It is worth noting that the complications related to (j2.74j) and (|2.75p arise only in the presence 
of several running coupling constants. Although beta function coefficients were studied for some 
models to a very large order, see e.g. [16], such computations typically involved only a single 
coupling constant. 



2.4 Perturbations by boundary changing operators 

Up to now we have assumed that there is a single (fundamental) boundary condition. The whole 
analysis can be easily generalised to the situation where we have superpositions of boundary condi- 
tions; in that case the set of boundary operators includes also boundary changing operators ipp b (x), 
where the two boundary conditions are labeled by a and b with a being the boundary condition to 
the left of x and b to the right. Local excitations of the pure boundary a are denoted The 
study of renormalisation group flows involving such operators was initiated in [23] . 

The OPEs of a bulk field approaching the boundary with label a, and that of two boundary 
fields have the form 



4>i(x + iy,x- iy) 



p "«Vr(lO + --- {y >x) . 



(2.76) 
(2.77) 



The only difference to the previous analysis is that there are now various superselection rules that 
demand for example, that products of boundary operators can only be non-zero if the intermediate 
boundary conditions match, or that the boundary fields that appear in the bulk to boundary OPE 
are always boundary preserving fields. Taking this into account, the boundary beta functions then 
have the following general form 



«y(«6) + a B P 5 ab x + pP(-V (aC V (Cfc) + Yl ^ V (a6) + • • • 



(2.78) 



where y^ = 1 — hp b are anomalous dimensions and are the coupling constants of the operators 

Tpf(x). The expression for the dimension shift (|2.16p generalises as 

anr btQr 
V p{aab) Z±_^ J)p(abb) U 



?p(ab) 



s 



p(ab) 



y, 



bb 



Similarly, the main results of the previous subsection (|2.7ip and (|2.72p now become 

L 

e jdx{<j>{x, e -)^ b {m h p a {oo)) 

-L 

' £\ Vr a 



S. 



p(ab) 



lim — 

e^O 2 



a id i 

D (aab)p D 4> 

' ' r P n.aa 
or 



pr 



b B r 

bb)p <P_ 

y bb 



?p{a.b) 



lim 

<5^0 



n—8 



^(<j6(e^X(0)^ a (oo)} 



ba i 



(2.79) 



(2.80) 



/ « °BI ( 1 

(aab)p <£_ I _J_ 

rp 2y™ \25 



^2 D (abb)p_l 



2yf \2S 



(2.81) 
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3 Some explicit examples 



Up to now our analysis has been very general. In this section we want to illustrate these general 
results with two simple examples. 



3.1 A single Neumann brane 

The simplest example is the case of a single Neumann brane on a circle of radius R. The action of 
this theory is simpl){^| 

(3.1) 



S=— cfzdXdX . 
2ttJ 

The bulk field that corresponds to changing the radius R is 4>(z,z) = 2dX(z)dX(z), which is an 
exactly marginal operator in the bulk. More specifically, we shall consider the perturbation 



5S = 2X dxdydX(w)dX(w) 



(3.2) 



that changes the radius R as R x = Re~ nX so that to the first order we have 5R = —ttRX. 

On the Neumann brane we have open string momentum states corresponding to the vertex 

5 and n € Z. We want to 



operators tp = e , whose conformal dimension is h = k with k = 
study how the conformal dimension of these operators changes as we change the radius. Thus we 
need to calculate^ 



£ = lim 

5^0 



2 r 5 d°d{e~ ikX (oo) dX(e M ) dX(e~ M ) e ikX (0) } - V] 



^z±_ ( _ 
25 



On the Neumann boundary we have dX = dX, and the correlation function equals 



2{e- lkX (oo)dX(e l »)5X(e 



1.0 \ 



-ii9\ JkX 



(0)} 



-2k 



1 



(z - zf 



(3.3) 



(3.4) 



Thus the integral is simply 

>*7V—5 

ikX 



d$ ( e~ lkX (oo) dX(e M ) BX(e- w ) e lkX (0) } 



JkX i 



-2 k TT - - cot 1? 
4 



-2k\ + — + 0{5) 
2d 



(3.5) 



The term that is singular in 5 is subtracted by the last term in (j3.3H . In fact, the only relevant or 
marginal boundary field that is switched on is the identity field with yo = I and = 1, and 

the corresponding bulk to boundary OPE coefficient is Bs 1 = 1 since 



2(dX(z)dX(z)) 



1 



1 

4y 2 • 



(z - zf 

Thus we find that £ = —2 k 2 it in this example, which implies that 

Sh = 2k 2 ir\ . 

This then agrees with the geometrical expectation since for h = k 2 with k = % we have 

Sh = -2k 26 -^ = 2k 2 7T\ . 
R 



(3.6) 



(3.7) 



(3i 



5 Throughout this section we set a' = 1. 



3 Note that tp is not a self-conjugate field, and we therefore have to insert the conjugate field at infinity. 
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3.2 Branes at angles 



A somewhat more interesting example is the configuration of two Dl-branes that stretch diagonally 
across a 2-torus, crossing each other at an angle (see figure 1). This brane configuration is obviously 
unstable since the relative open string between the two Dl-branes is tachyonic but this will not 
be important in the following. (One can imagine that this is only part of a more complicated 
background involving additional directions, and that the boundary conditions of these D-branes in 
the other directions are chosen so that the relative open string is not tachyonic.) 




Fig. 1: Radius perturbation on the torus: changing the radius of R\ modifies the relative angle 
between the two branes, and hence the conformal dimension of the corresponding boundary changing 
field. 



For simplicity consider the situation where the T 2 torus is orthogonal with radii R± and i?2- 
The bulk operator that changes either radius is an exactly marginal bulk operator, but it does 
have an important impact on the boundary theory since the ratio of the two radii determines the 
conformal dimension of the lowest string excitation between the two Dl-branes. In the following 
(section 3.2.2) we shall calculate the change in conformal dimension using the RG formalism we 
have developed above. As we shall see, this will reproduce the standard formula for the conformal 
dimension of boundary fields on branes at angles that will be reviewed in section 3.2.1. 



3.2.1 The geometrical analysis 

The torus theory is described by the action 

S = J (PzdX^dX^ . (3.9) 

We shall now consider changing the radius R\ by means of the perturbation 

5S = 2A jdxdydX 1 (w)dX 1 (w) , (3.10) 

which to first order gives 5R\ = —ttRiX. 

The two D-branes stretch diagonally across the torus; their angle relative to the x 2 axis will be 
denoted by ±0/2, where satisfies 

tan| = |i. (3.11) 
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The open string that stretches between the two branes satisfies the + boundary condition 

8X 1 {z) = -cosQdX l (z) + sin®dX 2 (z) , 

dX 2 (z) = sin e d X 1 (z) + cos edX 2 (z) (3.12) 
at one end (say for z = z on the positive real axis), and the — boundary condition 

dX 1 (z) = - cos6 ax 1 (z) - sine dX 2 (z), 

dX 2 (z) = -sinGdX 1 (z) + cosOdX 2 (z) (3.13) 

at the other (say for z = z on the negative real axis). By going to complex variables, i.e. by writing 
Z + = -^(X 1 + iX 2 ), Z~ = "^(X 1 — iX 2 ), we can write the open string fields as 

Z+(z, z) =i x Hy{ - a "~" e ~ ie , a ~ m \ v - ) (3.14) 

' mZ v" 1- "^"" (m + v)z m+u J v ; 

Z~{z,z)=ufly ( 7 a ™+" -j Q , °V ) , (3-15) 

where ^ = 0/7r € [0, 1). The modes a^i/ satisfy the canonical commutation relations 

[ a m-vi a n+p] = {™<- v)$m,-n , (3-16) 
and the Virasoro generators can be expressed in terms of them as 

L m = y : a+_ k _ u a~ +u : + \v{\ - v)5 mfl . (3.17) 

The conformal dimension of the lowest boundary changing operator ip h is thus 

h- + = - u) . (3.18) 

According to the analysis of [U [25] , the two D-branes will respond to the radius changing bulk 
perturbation (|3.1U|) by simply adjusting themselves infinitesimally, so that they continue to stretch 
diagonally across. To first order in A, the angle 6 thus changes via ()3.1ip as 

SQ = -7rsin0 A . (3.19) 

With ivv = O this implies that the conformal dimension of the lowest boundary changing operator 
changes as 

6hj, + = 2<>- l)sin9 A + C(A 2 ) . (3.20) 
This is the result we now want to reproduce using the RG approach explained above. 

3.2.2 The RG approach 

Formula (|2.8ip applied to the situation at hand reads 



£t { - +) = lim 



5^0 



5 



2 I d^(dX 1 (e rd )dX 1 (e-^)ip- + {0)ip + -(oo)) (3.21) 
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where the index r runs over all relevant boundary operators in the respective sectors. Since the 
model at hand is Gaussian the only relevant operator induced on the boundary by cj)(z, z) is the 
identity operator in the respective + or — sector. The corresponding bulk to boundary OPE 
coefficients can be read off from the expectation values 



2{dX 1 (z)dX\z))± 



<p_ 

4y 2 



COS 

4y2 



(3.22) 



which can be computed using the mode expansions (13.141) and (13.150 for v = 0. (The string fields 
in the presence of a single boundary are of the same form as (|3.14|) and (|3.15|) but with v = 0.) 
Thus 



cos . 



The three-point correlator in (|3.2ip is given by the one-point function of the radius changing op- 
erator in the presence of the boundary conditions (|3.12p and (|3.13p . A straightforward computation 
yields 



2(dX 1 (z)dX l (z)tl;- + (0)^ + -(oc)) 



2(dX 1 (z)BX 1 (z)) e 
1 _ ie z» z(l-v) + zv 



+ c.c. 



z[z — z 



(3.24) 



Viewed as a function on C rather than on H + , the correlator has a logarithmic branch cut along 
the negative real axis. The integral at hand can be easily evaluated 



TV— 8 



2 J di9{dX 1 (z)dX 1 (z)) 



e 



p -iB+2ivd 
d& - . „ : (1 - v + ve-™) + c.c. 



8sin 2 tf 



(3.25) 



cos(-G + <&(2u - 1)) 



4sin"!9 



l9 = 7T-5 



■9=8 



cos(-e + ^(2t/ - 1)) 
2 sin 5 



where in the last step we used = irv. Substituting (|3.25|) and (13.230 into (|3.2ip we obtain 



£. 



4i~+) 
<f>4> 



lim 

8^0 



cos(-9 + 5(2u - 1)) cosG 



2 sin 5 



+ 



Noting that 

6h ~4> + = - S y^~ 

we finally get the same result as in (I3.20p . 



-\£ 



26 



(2v - 1) sinG 



(3.26) 
(3.27) 



4 Third order coefficients in the pure bulk or boundary case 

In section 2 we explained how to calculate higher order coefficients in a theory with bulk and 
boundary perturbations. Actually, the techniques used there can also be easily generalised to the 
pure bulk or pure boundary case; this will be sketched in the following. We begin with a discussion 
about which coefficients of the third order terms in the beta functions contain universal quantities, 
paralleling the discussion in section 2.1. In section 4.2 we then describe how to obtain useful 
formulae for these coefficients. 
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4.1 Universal quantities 

In this subsection we shall only consider the pure bulk theory; the discussion for the pure boundary 
case is very similar. Consider a conformal field theory perturbed by 



5S = Y^Xl~ Vi jd 2 z (j)i{z,z) , (4.1) 



where as before yi = 2 — Ai are the anomalous dimensions, I is a renormalisation length scale, and 
Aj are the dimensionless coupling constants. The beta functions have the general form 

P l = Vi% I Y. Cl U XXI ■ E4* A<A ' A * + 0(A 4 ) . (4.2) 

ij ijk 

We take the constants C\j and J 7 ^ to be totally symmetric in i,j and i,j, k, respectively. Under a 
general change of scheme the coupling constants are redefined as 

A' := A' + £ <V,A'A' + £ // ;/ ,A' A ; A A ' + 0(A 4 ) , (4.3) 

ij ijk 

where the ■ are, without loss of generality, symmetric in i and j. The beta functions in the new 
scheme are 

$ = yi \ l +Y / ~X i ~X j (cl j + 4 j (y i + y j -y l )) 

+52%% \" P tjk + f l ijk ( yi + Vj +y k - m) (4.4) 

ijk 

+ 3 Em Sperm(ij, fc) ( C mfijk ~ ^mi C jk ~ C mi C jk(y m V* ~ Vl) 



+ (9(A 4 ) . 

We observe that the second order coefficients C\j do not change under this transformation if and 
only if the second order resonance condition yi + yj = yi is satisfied. As for the coefficients F\j k 
at the cubic powers of the couplings, it can be seen from (|4.4p that the basic requirement for F\- k 
to be universal is that it satisfies the resonance condition yi + yj + y k = y\. However, even if the 
resonance condition is satisfied, the third line in (|4.4jl shows that the corresponding coefficient may 
not be invariant under general scheme changes because of the lower order coefficients C\y The 
resulting transformations of the resonant coefficients are parametrised by the tensors c£-. For an 
n-dimensional coupling space the dimension of the space of coefficients J~\j k is n ( w +i)( n + 2 ) w hii e 

that of the coefficients c^- is HiB+li_ Depending on how many coefficients are resonant, there may 
be some functions defined on these resonant coefficients which are invariant and thus give universal 
quantities. For example if all couplings are marginal (y^ = for all i) then generically there must 
be a subspace of scheme independent coefficients of dimension n < - n+1 ^ n ^ m 

While in general it is hard to write out explicit expressions for universal quantities in terms of 
T\ai. and £■„■ we can do so in the absence of second order resonances because we can then use a 
special scheme in which all C\j vanish. Given a cubic resonance yi + yj + y k = yi, the values of the 
cubic coefficients F\j k in that scheme are universal and can be expressed via the coefficients in an 
arbitrary scheme as 

f nk = f ijk + o ■ • (4-5) 

perm(ij,fc) «i 

The scheme independence of (|4.5|) can be checked directly using (|4.3|) and (|4.4p . 
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We can also consider a situation analogous to the one considered in section 2.1 when the universal 
quantity gives dimension shifts under a truly marginal deformation. Let A be a coupling constant 
corresponding to an exactly marginal operator <p(z,z). The beta functions for the other operators 
(f>l(z, z) have the form 

k ij 

+ £ r'^xx'x 1 - + \ 2 e l H + \ 3 j%+ + . . . . (4.6) 

ijk 

We can make the beta functions ft 1 homogeneous in X 1 up to cubic order by a coupling constant 
redefinition 

ji = X i + C Jtt\ 2 + 5^ A 3 . (4.7) 

Vi Vi 

This redefinition is possible because A is truly marginal. The anomalous dimensions of the operators 
<f>i are then given by the eigenvalues of the matrix 

DlW - [ ?| | (4.8) 



\ k =o 



and have the form 

y i [X]=y i + \8 l p + \Hf ] + ... . (4.9) 
A straightforward computation yields 

<f } = c;, (4.io) 

and 

pk pi pk pi 

f (2) _ -r-i ^tptj^ik sr^ ^<j>&<t>k 



x\ 2 ) _ -r-i o \ " 99 tk , \ " <P» fc / yi -i i \ 



Vk rr. Vi - Vk 

k k^pi 

The coefficients are resonant and thus universal. One can also check that (|4.11j) is invariant 
under an arbitrary coupling constants redefinition of the form 

\ l = f\\) + ^fl(\)\ k + ^fL(\)X\ k + ^rt jk \ l V\ k , (4.12) 

k ik ijk 

where / , f l k , f\ k are polynomial functions of A. 

Finally let us mention the well known fact that if there is a single running coupling constant 
whose UV dimension is marginal, then both the quadratic and cubic terms in its beta function are 
universal. 



4.2 Computation of coefficients 

Now that we have understood which coefficients are universal, we can ask how they can be calculated 
explicitly. As in the bulk-boundary case discussed in section 2, we can either use a minimal 
subtraction scheme (see section 2.2) or the OPE scheme of section 2.3. As before, the minimal 
subtraction scheme is conceptually clearer since one does not need to introduce beta functions for 
irrelevant fields. However, the calculation is somewhat unwieldy in this scheme, since one has to 
isolate the divergences in the UV cut-off e for finite IR cut-off L. 

In the following we shall only consider universal quantities for which the calculation in either 
scheme must give the same answer. Since the OPE scheme is technically simpler, we shall use it 
to determine explicit expressions for these coefficients. We have also checked that our result agrees 
with what would have been obtained in the minimal subtraction scheme (as must be the case). 
Moreover for brevity we will focus on the quantity (|4.5H . It is straightforward to extend our results 
to the dimension shifts (|4.11|) and to a cubic term in a beta function of a single marginal coupling. 
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4.2.1 Resonant bulk coefficients 



In the OPE scheme the RG equations are determined from the condition that the variation ed € e ss 
vanishes in the limit e — > 0. As before we regularise the theory by point splitting, i.e. we introduce 
a sharp UV cut-off e. In addition we have an IR cut-off L. To cubic order in the couplings we have 



1 + ]T \ i e~ m fd 2 z Mz, z) + 1 Y. X %< ''' "' j [d 2 zid 2 z 2 9 12 <j> i (z 1 ,z 1 )<l> j (z 2 , 

i ' ij 

+^J2 Xixjxke ~ Vi ~ Vj ~ m [ [ [d 2 z l d 2 z 2 d 2 z 3 9 u 9 23 e 1 ^ l (z 1 ,z 1 )<P J (z 2 ,z 2 )4> k (z 3 ,z 3 ) + ... , 



ijk 



where i 



= 9(\zi — Zj\ — e)9(L—\zi — Zj\). The variation ed e of this expression can be computed using 
Setting ed e e ss ~ at second order in the couplings one obtains the well known expression 



nm (~* m 

ij — ij i 



(4.13) 



where C 4 - are the bulk OPE coefficients (j2.2|) . At the cubic order we have the equation 



~ A*A J 'A* V 

perm(i ,j,k) 



1 



e -w-w-»*+i / Sz x d 2 z 2 d z z 3 5\ 2 9 13 9 23 ^{z x ,zMz 2 , z 2 )cj> k (z 3 , z 3 ) 



d zid z 2 9 12 <p m (z 1 , Zl)(j) k (z 2 , z 2 ) 



(4.14) 

As was discussed in section 2.3 the above equation in general may still have divergences. However, 
in the resonant case, i.e. if yi + yj + y k = yi no such complications arise. Then we can write 

(Hjklres = Ymi-^ Jd 2 z 2 9 20 ^ J d 2 z 1 9 w 5l 2 {<pi(zi,z 1 )(p j (z 2 ,z 2 )(j) k (0)<j)i(oo)) 



perm(ijk) 



(4.15) 



where 9±q = 9{\z\ \ — e)9(L — \zi\), and similarly for <5io- As before, we consider spinless fields, for 
which we can express the four-point correlator in terms of conformal blocks, 



i(zi,Zi)fa(z>i,Z2)<f>a(z3,Z3)fa(z4,Z4)) = Yl\zij\ 2(S hi hj) Y 12U (r],r)) 



i<j 



with 



and 



1 4 

Z 3 ' 6= oY. Hi ' 



i=l 



Zl 2 Z 3 4 
Zl 3 Z 2 4 



{ J: m Ci2 m c m3 *F% M (n)F% M (fj) 

£ m C 32 m C ml *F£ M (l - V )FE M {1 - fj) 

k E m Ci3 m C m2 4 F^ 2 ,(l/v)Fr 3 M 1 /^ ■ 
The conformal blocks are normalised such that 

*i2,34fa) ~ v hm ~ 5 Ko 3 M ~ fj h ™- 6 for n, fj -> o. 



(4.16) 



(4.17) 



(4.18) 



(4.19) 
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Replacing the variables in the first part of the square brackets in (|4.15p by an angular variable ip 
and the cross ratio 77, 

77 

the angular variable can be integrated out by the integral over z±, and one finds 
(4- fc ) res = }mljd 2 v £ e(l-\ V \)e(\ V \-i)U±-Rer ] )6(\i 1 \ 2 (£-l) + 2Re V -l) 

perm(i,j,k) 



TO >' 



In the second term in the bracket of (|4.15l) . we changed variables to Z2 = e/??. The function in the 
last line in (|4.20p can be integrated explicitly, and we obtain 



(Fijk)res = Jim 

perm(i ,j,k) 



d 2 r,9(l-\ V \)e(\rj\-i)e(l-Re V ) 

7T 2 ^ Cij m C m k / e\Vl-Vh-Vm 



W -*-«-ai*-4A44| 1 _ vl 2A + y J+ y k -4 Y ■ { j _ « £ 

3 m m ~ yk ~ ym L 
_2 ri mri I 

+y £ ^ ^ . (4.21) 

The universal quantity defined in (|4.5p is then simply given by the first two lines of (|4,2ip . 

It can be checked using the asymptotics (14. 191) and the properties of conformal blocks (|4.18ft that 
the integral (|4.2U|) converges in the regions 77 ~ e/L — > 0, |1 — 77 1 ~ e/L — * and \r/\ ~ L/e — ► 00. 
One can thus safely set e = in (I4.20p to obtain an integral expression 



perm(i,j',fc) 

-£)Cii m C TO fcV i ^~ 1 ' m ~ 2 } • (4.22) 
m 

Note also that expression (|4.22p is L independent. In particular this means that the infrared 
divergences that were present in individual summands in (I4.21|) mutually cancel each other. This 
agrees with the general results of [T71 [T8J 

By suitable changes of the integration variable 77 in the terms with permuted indices k 
it is possible to write fL^ by means of integrals over three disjoint subsets tiling the whole 77- 
plane. Consider the transformation 77 1— ► 1 — 77, for which the cut-off functions in the integral (|4.22|) 
become 9(1 — \rj — l|)0(Re?7 — 5). The asymptotics (|4.18p for Yij^i are such that the divergence 
of the transformed integrand that arise from the limit 77 1 — >- 1 is again canceled, once we take the 
transformed subtractions, i.e. the second line in (|4.22[) . into account. The other transformation is 
77 1— > I/77. In this case the cut-off functions read 9(\rj\ — l)9(\rj — 1| — 1) after the transformation, 
and the divergence of the corresponding integrand for 77 — > 00 is canceled as well. Together, the 



7 The perturbation expansion for Wilson coefficients proposed in [171 118] was shown to be IR finite to all orders 
under certain assumptions on the UV renormalisation scheme. As we are interested in scheme independent quantities 
their result applies. 
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regions carved out by the cut-off functions for the three coordinate choices tile the whole rj plane. 
Using this we can recast (|4.2'2|) as 



[HjkU = ~ yrfN[l^| 2r+K+% - 4 |l-r/| 2r+%+ ^-%-,H(r/,f/)-5^(r ? )J , (4.23) 

SuM = ^c ij m c mk l \ r] \y* + y^- 3 e(i-\ri\)e^--Rer,) 

m 

+ Y c kJ m c mi l |i - r)\y« + yi-y™- 2 e{i -\ v - lQ^Rer? - §) 

m 

+ ^C lk m C m / M-" +2/m ~ 2 0(M - 1)0(\V - 1| - 1) • (4-24) 

m 

In this form the integration runs over the whole 77-plane. Although the subtraction function Sijki(jj) 
still has a piecewise form it is expressed quite explicitly. 

4.2.2 Resonant boundary coefficients 

On the boundary the computation can be done in a similar way as in the bulk. We consider a 
boundary perturbation of the form 

5S = Y M* t~ Vs [ dx i) s (x) , (4.25) 

where now y s = 1 — h s . Up to the third order in the couplings the RG equations take the form 

(X s = y s ^ s + Y ?>>V + £ S&r^W + • • • • (4-26) 

p,q p,q,r 

As before, we only introduce counterterms at the quadratic order for marginal or relevant fields, 
and the corresponding coefficients are 

T>U = D pq s , (4.27) 

where D pq s is the OPE coefficient of two boundary fields (|2.3p . In the resonant case where we have 
Vs = Up + Vq + Ur, the coefficient (Gp qr )res can be written as 

(0 s pqr U = ^lim Y { e J 2 (^(O)VffWr(^(oo)) 



where 



perm(p,g,r) 



L (4.28) 



_ ^ y j ^ p ( a; )^(o)^(oo)) I . 

By similar arguments as in the previous subsection we can find 

\J dr > E {v r+yp+yq - 2 a-vY +Vq+yr - 2 Yp q ,rs(v) 

Y DpfDtr'rpm-vt-i - Y WC 1 - ^+^-1 j 



(^p(jr)res 



perm(p,<7,r) 



(4.29) 
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where 

y pq ,rs(v) = £ D P« tD tr S Fp q ,rs(ri) = £ D pt s F^ (1 - r?) . (4.30) 



Here the conformal blocks F p {rj) have cuts running from — oo to zero, and from 1 to +oo. In 
addition, their asymptotic behaviour is 

F^V) ~ r] h *- s (4.31) 

where 5 is defined as before, i.e. S = \{h p + h q + h r + h s ) . Finally, the scheme-independent quantity 
is given by 

Q s -(Q s ) +- V ST D P<i t ( D tr S + Drt s ) ( , 

y pqr -{y pqr )™ + 6 L y s -y r -y t " ^ 

perm(p,q,r) t 

In the case where several irreducible boundary conditions are involved, one has to keep track 
of their labels, and bear in mind the superselection rules, in particular the order of operators. 
This leads to additional splittings and recombinations of the integrals over four-point functions and 
subtractions. Apart from this technicality, it is however straightforward to include the boundary 
labels. We have refrained from writing them explicitly to keep the formulae simpler. 



5 Conclusions 

In this paper we have studied conformal perturbation theory beyond the leading order. We have 
shown that, at least up to quadratic order, the combined bulk boundary perturbation problem is 
renormalisable, using the minimal subtraction scheme. We also discussed the more commonly used 
'Wilsonian' OPE scheme, and found it to have some shortcomings at higher order in perturbation 
theory. We identified systematically the universal (scheme-independent) quantities, and gave ex- 
plicit formulae for them at third order in terms of integrals of conformal 4-point functions. Finally, 
we explained how essentially the same analysis works for the pure bulk and pure boundary case. 
It seems plausible that similar techniques should allow one to prove renormalisability at arbitrary 
order in perturbation theory, but we have not attempted to do so. 

Our work was originally motivated by the question of how the dependence of the conformal 
dimension of a boundary changing field upon a bulk modulus can be understood from the world- 
sheet perspective. Our considerations demonstrate that this effect is captured by a certain universal 
quadratic RG coefficient, for which we gave an explicit formula. This result should also have 
interesting applications in other contexts; in particular, it provides a world-sheet method to study 
the stability of brane setups under arbitrary bulk deformations. 
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